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Abstract 

Trajectories of charged particle in combined poloidal, toroidal magnetic field and 
rotation-induced unipolar electric field superposed in Schwarzschild background 
geometry have been investigated extensively in the context of accreting black 
holes. The main purpose of the paper is to obtain a reasonably well insight on 
the effect of spacetime curvature to the electromagnetic field surrounding black 
holes. The coupled equations of motion have been solved numerically and the 
results have been compared with that for flat spacetime. It is found that the 
toroidal magnetic field dominates the induced electric field in determining the 
motion of charged particles in curved spacetime. The combined electromagnetic 
field repels a charged particle from the vicinity of a compact massive object and 
deconfines the particle from its orbit. In the absence of toroidal magnetic field the 
particle is trapped in a closed orbit. The major role of gravitation is to reduce the 
radius of gyration significantly while the electric field provides an additional force 
perpendicular to the circular orbit. Although the effect of inertial frame dragging 
and the effect of magnetospheric plasma have been neglected, the results provide 
a reasonably well qualitative picture of the important role played by gravitation 
in modifying the electromagnetic field near accreting black holes and hence the 
results have potentially important implications on the dynamics of the fluid and 
the radiation spectrum associated with accreting black holes. 
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1. INTRODUCTION 

From the nearest stellar object Sun to the furthest cosmological object quasar, the pres- 
ence of the magnetic field effectively governs not only various physical phenomenon but also 
the evolution of many astrophysical objects. 

The influence of the magnetic field on the dynamical flow of the infalling plasma around 
a compact object and on the emergent radiation spectrum have long been realized in the the- 
oretical study of very high energy astrophysical objects like pulsars, quasars, Active Galactic 
Nuclei (AGN) and X-ray binaries. Unfortunately no convincing theory has been well estab- 
lished so far which can explain the radiation from these sources consistently and correctly. 
Though the definite mechanism has not been understood properly, it is generally believed 
that the radiation emission could be due to the plasma processes near compact objects like 
neutron stars and black holes. Thus it becomes necessary to consider plasma processes in 
intense-gravitational-field backgrounds, as such compact objects would be massive enough to 
produce significant spacetime curvature effects. Therefore one has to consider the dynamics 
of the fluid in curved spacetime which is quite an involved problem since one requires to 
solve the general relativistic magnetohydrodynamics equations. 

It is well known that the best way to understand the structure of any field is to study 
the dynamics of test particles in that field. If a particle is charged then it deviates from 
its geodesic motion and the study of such trajectories would reveal informations about the 
influence of the interacting field on the spacetime geometry and vice-versa. The most im- 
portant in the astrophysical context is the electromagnetic field. In an electromagnetic field 
a charged particle gyrates with a definite radius of gyration and that by and large deter- 
mines the radiation spectrum of the object. The concept of radius of gyration in classical 
mechanics is modified when one considers curved spacetime 1 . Hence it is very important 
to investigate the trajectory of charged particles in curved spacetime with different types of 
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electromagnetic field configuration. 

For this purpose one should look for solutions of Einstein-Maxwell equations which are 
asymptotically flat and have non-zero dipole magnetic field even in the absence of rota- 
tion. Although these systems of equations are formidable to solve in general, there are some 
solutions obtained by perturbation technique under the assumption that the electromag- 
netic field is weak compared to the gravitational field and thus it does not affect the basic 
geometry. This is achieved essentially by solving the Maxwell equations on a given curved 
spacetime background geometry. This is indeed a valid assumption for even the most intense 
magnetic field associated with pulsars carries an energy which is very small compared to the 
gravitational potential energy on the surface of a neutron star of one solar mass. Therefore 
it is considered that the magnetic field would not affect the spacetime curvature but the 
curvature could affect the magnetic field. 

With such a consideration Ginzburg and Ozernoi 2 , Petterson 3 , Bicak and Dvorak 4 have 
obtained solutions for dipole magnetic field on Schwarzschild background, whereas using a 
similar approach Chitre and Vishveshwara 5 , Petterson 6 have obtained solutions for stationary 
electromagnetic fields on Kerr background. 

Charged particle dynamics in such electromagnetic fields have been extensively studied 
by Prasanna and Varma 7 for the Schwarzschild background, by Prasanna and Vishveshwara 8 
and by Prasanna and Chakraborty 9 for the Kerr background. 

Prasanna and Sengupta 10 have obtained trajectories of single charged particle in the 
presence of toroidal magnetic field superposed on the Schwarzschild background geometry. 
In that work it is shown that the toroidal component of the magnetic field repels the particles 
away and particularly the incoming ones get ejected in jet-like trajectories. Hence it has been 
pointed out that to produce jet like features from the accreting matter in the high energy 
astrophysical scenarios, it is very essential to have the toroidal component of the magnetic 
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field associated with the central gravitating body. 

It is known that almost all of the celestial bodies do have a non-zero angular momentum 
and therefore one must consider the important role played by the angular velocity of the 
central object in determining not only the spacetime geometry exterior to it but also the 
electromagnetic field configurations. For this purpose Kerr metric has been used in the 
earlier works and the motion of the particle in the equatorial plane has been investigated. 
Since the exterior region in these studies is considered to be absolutely vacuum, the electric 
field induced by rotation of the object is obtained to be quadrapolar in nature. However this 
cannot be a realistic situation as it is well known that in the case of pulsar magnetosphere 
and in the accretion disc the region in which the motion of a charged particle is considered, 
are plasma filled. As a consequence Ohm's law should be taken into consideration within the 
entire region while determining the electric field. Due to the Ohmic induction a rotation- 
induced unipolar electric field must be present in that region if one neglects the higher 
multipole by considering comparatively slow rotation. 

In the present paper I have presented the trajectories of a charged particle under the 
combined poloidal, toroidal magnetic field and the rotation-induced unipolar electric field 
superposed in Schwarzschild background geometry. This is an extension of the work by 
Prasanna and Sengupta 10 . The main purpose of the paper is to obtain a qualitative picture 
of the role played by toroidal magnetic field in curved spacetime in the vicinity black holes. 

2. THE ELECTROMAGNETIC FIELD IN CURVED SPACETIME 

2.1. The spacetime metric 

Since we are dealing with the exterior geometry of black holes we need the vacuum 
solution to Einstein's field equations as the self gravitational field of the surrounding plasma 
is negligible. 

For a rotating compact object one may consider the Kerr metric which has the following 
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form 



ds 2 = 




c 2 dt 2 - 2a 



2 2mr sin 2 9 



cdtd(f) + 




-a 2 A sin 2 9 



sin 2 9d(f) 2 + 



'±dr 2 + p 2 d9 2 , 



(1) 



where 



A = r 2 — 2mr + a 2 , p 2 = r 2 + a 2 cos 2 9, a = J/m, 



J being the total angular momentum and m = -Jr, where M is the total gravitational mass 
of the object. 

For the case of a black hole the effect due to the rapid rotation upon the spacetime 
curvature could be important and the Kerr metric may be more suitable. Nevertheless if one 
considers black holes with ^ >> j| one can adopt the Schwarzschild metric, for small j| the 
Kerr metric reduces to the Schwarzschild metric and the problem becomes much easier to 
tackle with under this approximation. In the present investigation we consider sufficiently 
slow rotation of the object so that the effect of rotation on the background geometry is 
insignificant. Hence we consider Schwarzschild background geometry. The investigation 
with Kerr background geometry is under progress. 

The Schwarzschild metric is given by : 



where m = ^r, M being the total gravitational mass of the object. 
2.2. The magnetic field 

The vector potential for the poloidal magnetic field in Schwarzschild background geometry 
is given by Petterson 3 , Prasanna and Varma 7 , Wasserman and Shapiro 11 as 



(2) 





[r 2 ln(l ) + 2mr + 2m 2 ]. 



(3) 



8m 3 
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Using the definition = (Aj ti — A it j) we get the components of the magnetic field in 
Schwarzschild background geometry as 

„ 3asin 2 # r r , . 2m, , 2m, , , 

F<t>r = -. 2 ~ In 1 + 1 ) +1 , 4 

4m 2 m r r 

„ 3ur 2 sin 9 cos 9 r , . 2m. 2m. m xl 

^ = , 3 In 1 + — 1 + - , 5 

4m d r r r 

being the dipole moment associated with the central source. 
In magnetized accretion disks around a compact object the surrounding plasma is sup- 
posed to be tied to the magnetic field and its inertia would cause the magnetic field lines to be 
bent backward creating a toroidal component. The toroidal field component in Schwarzschild 
background geometry is given as 10 

Fre = r£ (r - 2m) ^ 
r D (r — 2m) sin# 

r D and B being constants. 

2.3. The induced electric field 

We consider the electromagnetic field to be stationary. For this we require the magnetic 
dipole axis to be aligned with the rotational axis of the central object. In the present study we 
ignore any effects due to the surrounding plasma. Since the surrounding medium exterior to 
the object is considered to be plasma filled, there will be a rot at ion- induced unipolar electric 
field such that E.B = everywhere. The components of the induced electric field considered 
here are obtained by using the generalized Ohm's law (assuming infinite conductivity of the 
medium) 

i^V = o, 

where 0, 0, v?) are the components of the four velocity vector of the corotating 

plasma. 



In Schwarzschild background geometry the components of the induced unipolar electric 
field are given as 12 

_, 3uu sin 2 9 r 2m. 2m. , 

= — - Mi + i + 1 , 7 

4m^c m r r 

„ 3uur 2 „ „ r , , 2m. 2m. m xl . . 

= -r^-sin0cos0pn(l +— 1 + - , 8 

4m d c r r r 

where uj is the angular velocity of the central object as measured by an observer at infinity. 
From equations (7) and (8) we obtain the potential for the electric field as 

3uur 2 sin 2 9 r , , 2m. 2m. m xl 
A = -^— 3 In 1 + — 1 + - ■ 9 

The electric field arises due to the corotation of the plasma with the object and is con- 
tributed by the poloidal magnetic field. The toroidal magnetic field has no contribution to 
the electric field as can be seen from the generalized Ohm's law. 

Therefore in the present discussion we are considering a region with combined poloidal 
and toroidal magnetic field and rotation-induced unipolar electric field superposed on 
Schwarzschild background geometry. 

3. THE EQUATIONS OF MOTION OF A CHARGED PARTICLE 

3.1. In Schwarzschild background geometry 

For the motion of a charged particle of charge e and rest mass m Q there are two constants 
of motion, the canonical angular momentum / and the energy E as given by 

V<t> + eA^ = /, (10) 

v t + eA t = -E, (11) 

wherein all quantities are normalized with respect to the particle rest mass m Q . The velocity 
four vector of the charge particle, v l is timelike, which for the metric of signature +2 is 



expressed as 



g ijV % v 3 = -1. 

From equation (10) we get 



dd> I 3eu ri ,„ 2m. 2m ,„ m x . 



From equation (11) we get 

. cdt ,„ 2m. 3euu;r 2 sin 2 # ri . 2m. 2m ,„ m N1 , 
u* = — = (1 ) {E + p 5 In 1 + 1 + — ) }. 13 

From the covarient Lorentz equations 

v!jV j = eFjV j , 

where v l — ^ is the four velocity of the particle and the semicolon denotes the covariant 
derivative with respect to the spacetime metric, we obtain 

d 2 r m. 2m. i.dr, 9 . 2m,,.d# x9 9 n .dd). ^ 

ra,„ 2m..cdt. 9 e 2m. , rB (r — 2m) d# 

— (1 )( — ) H (1 ){ 

r 2 r rfs m D c 2 r r D (r — 2m) sin 6* ds 

3/iSin 2 2 V --) + (—-!) ln(l - ^)](1 - _ (14) 

2m 2 LV r V 2m 7 V r m r ; K ds ds Jf ' V 7 

d 2 # 2drd6 . n n/ d(j), 2 e r rB (r -2m) dr 

— — ~ = — - + sin0cos0(— -) — — - — r —- 

ds 2 r ds ds ds m c 2 r 2 r (r — 2m) sin 9 ds 

3usin# cos# , 9l .„ 2m. „ . dt dd>., . . 

4m' {r ln(1 " - ) + 2m<r + m)}( ^ " dl' 1 ' (15) 

3.2. In flat spacetime 

In order to investigate the effects of spacetime curvature we need to determine the tra- 
jectories of the charge particle in flat spacetime as well. These will provide the motion of a 
particle in the Newtonian mechanics under various electromagnetic conditions. However in 
the present discussion we neglect the toroidal magnetic field component for flat spacetime 
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since the poloidal magnetic and the induced electric fields are sufficient for a comparative 
study of the particle motion in curved and flat spacetimes. 

Expanding in Taylor series and neglecting the higher order terms containing m we obtain 
from equations (12) to (15) the equations of motion of a charge particle in electromagnetic 
field without gravitation as 

d 2 r f( d9 2 . 2Q J(f) 2 efism 2 6 d<f> dt 

- = r{(-) + sm *(-) } - — ^ - u-) , (16) 

d 2 9 IdrdQ . Q(^) 2 2e/isin#cos# dt d<p 

ds 2 r ds ds ds m c 2 r 3 ds ds ' 

d(p I e/i . . 



ds r 2 sin 2 6 m c 2 r 3 ' 
cdt efiusin 2 9 

— = k, — . (19) 

ds m c 6 r 

4. NUMERICAL SOLUTIONS 

As the orbit equations involve transcendental functions it is analytically impossible to 
make any general analysis. Hence we restore to numerical integration and get the quali- 
tative picture of the nature of trajectories for different values of the physical parameters 
appearing in the equations. For convenience of calculations we shall consider the equations 
in dimensionless form by introducing the dimensionless quantities 

r s I c t muJ x e B Q eu 
p= — , a = — , L = — , r = — , W = , A T = -, A P = — — . 

m m m m c m c z m c A m A 

With this definition the equations of motion for curved spacetime read as 

-In - h(-v + AAp(1 " ^ dd _ 

p 21 p ,[ da ) sm9 da 

(20) 
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da 2 pdada da p 2 (l — ^) sin 9 da 

3Ap . „ ,„ 2. 2_ L, Iir dT dd). . . 

+— -sin0cos0{ln 1- - ) + - 1 + - } (W- f-) , 21 

4 p p p da da 

d<P L 3A Pri „ 2 X 2. 1,. 
acr /sinfl 8 p p p 

1 _ )-i {E + sin 2 e ln(1 _ + x + } (23 ) 

da p 8 p p p 

where A = 

Similarly, the equations of motion for flat spacetime read as (without the toroidal field 
components) 

= P{(f- f + sin 2 eA* } - % sin 2 9{f - W^) , (24) 
da 2 da da p 2 da da 

d 2 8 2 dpd9 . . n ,d(b,„ 2\ P . n n/TTT dr deb. . . 

+ sin#cos#(--) 2 f sin#cos#(W— - - --) , (25) 



da 2 p da da da p 3 da da ' 

ddj L \p 

da p 2 sin 2 9 p 3 ' 

dr „ X P W . 2 



(26) 



— = E — ^— sin 2 9. (27) 
da p 

In order to solve the equations of motion we need to pescribe the initial conditions. 
We assume initially p = p Q , 9 = 9 Q , <fi = <p and r = t d . Also at p Q , ^ = (^)o, % = 

\dcr)°- 

Throughout our discussion we take (^)o = and r = 0. Therefore, from equation (22) 
we get 

L = sin 2 (£ )[ln(l - A) + ^-{l + H (28) 

O Po Po Po 

From equation (23) we get 

(£). = (1 - + 5^ S ^(^)[ln(l - 1) + 1(1 + I)]}. (29) 

dcr p 8 p po Po 
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From the spacetime metric (2) we get 

<£>- = ± 00& ; -'< ): + 1}1I/ '- < 3 °> 

Equations (28) - (30) provide the initial conditions for the motion of the particle in curved 
spacetime. 

Similarly, from equation (26) we obtain 

L = —sm 2 (6 ). (31) 

Po 



From equation (27) we get 



A =B _^ sinW (3 2 ) 
da p 

From the Minkowskian metric 

d g> = _ c 2 dt 2 + dr 2 + r 2 dQ 2 + r 2 ^2 ^ 

we get 

(g )o = ± [i + (*);- rf (*)r. 03) 

Equations (31) - (33) give the initial conditions for the motion of the particle in flat 
spacetime. 

The equations are numerically solved by using the 4th order adaptive step size Runga- 
Kutta method. The integrations are performed between finite limits with an accuracy of 
1CT 7 — 1CT 8 , while changing the step size appropriately through self inspection depending 
upon the initial conditions. The trajectories are obtained in spherical polar co-ordinates 
which are then converted into Cartesian co-ordinates by using the relations 

x = r sin 9 cos 0, y = r sin 9 sin 0, z = r cos 9. 

The results are presented graphically. 
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5. RESULTS AND DISCUSSIONS 

As obtained by Prasanna and Varma 7 , when the magnetic field is large, the motion of a 
charged particle looks similar to that in the flat geometry wherein the particle gyrates in a 
given tube of lines reflecting between two mirror points located symmetrically with respect 
to the equatorial plane (see Figure la). On the other hand, if the magnetic field is lower, then 
the particle oscillations up and down the equatorial plane damp continuosly as r decreases 
and eventually the particle gets sucked in by the central object. The inclusion of the induced 
unipolar electric field changes the nature of the orbit by introducing a E x B/£? 2 drift. The 
particle still remains trapped in a mirror configuration reflecting between two points on 
either side of the equatorial plane (see Figure lb) just like the case without the electric field. 
Figure lc shows the trajectory of the particle in a combined poloidal and unipolar electric 
field superposed in flat spacetime. A comparison between Figure lb and lc reveals that 
the gravitational field mainly contributes in decreasing the radius of gyration of the particle 
significantly i.e., in curved spacetime the particle gyrates in much tighter orbit than it does 
in flat spacetime. 

From equation (12) to equation (15) we notice that unlike the case of pure poloidal field, 
a particle on the plane 6 = tt/2 having zero initial velocity in the meridonial direction is 
not confined to the equatorial plane as the acceleration along the 9 direction is still non-zero 
in the presence of toroidal magnetic field. The inclusion of rotational motion of the object 
modifies the initial conditions. From the initial conditions we notice that the initial velocity 
of the particle is well defined even if E = since the particle acquires its initial velocity 
due to its interaction with the electric field. The situation is similar to the case for the Kerr 
background geometry 8 . But here we have considered a rigid rotation neglecting the dragging 
of inertial frame effect. Since it is difficult to obtain the trajectory of a charged particle in the 
off equatorial plane with the Kerr geometry due to the presence of non-zero off diagonal terms 
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in the metric, the particle trajectory has not been investigated in that situation although 
the situation is more realistic in that case for astrophysical interest. The present study can 
very well provide a nearly similar qualitative picture with two differences i.e., the absence of 
the dragging of inertial frame effect and the rot at ion- induced quadrapolar electric field. As 
mentioned in the introduction, if we consider comparatively slow rotation, the inertial frame 
dragging becomes negligible and the consideration of rot at ion- induced unipolar electric field 
is much more realistic in the astrophysical situation than that of qudrapolar electric field. 

Prasanna and Sengupta 10 investigated the trajectories of a charged particle in combined 
poloidal and toroidal magnetic field superposed in Schwarzschild background geometry. The 
effect of rotation and hence the electric field has not been considered in that study. In that 
work it is found that when the toroidal field is extremely small compared to the poloidal 
field the particle trajectory appears somewhat similar to that in the case of pure poloidal 
field till the particle reaches a distance two Schwarzschild radii away from the central object 
where gravity seems to dominate and pull the particle in. As the magnitude of the toroidal 
field is increased, the nature of the trajectory alters and when the magnitude of the toroidal 
field becomes significant compared to that of the poloidal field the particle gets deconfined 
even very near to the central gravitating object and gets pushed away to infinity. The most 
interesting feature is that if the initial velocity of the particle is towards the central object 
then in the presence of the toroidal field the particle comes very near to the object and then 
gets ejected away to infinity. 

As the earlier investigation has not included the electric field induced by rotation one 
may expect significant changes in the trajectories in more realistic case. But here we show 
that the above situation remains the same even if one includes the electric field. Figure 
2 and Figure 3a show the trajectory of a charged particle in combined toroidal, poloidal 
magnetic field and induced unipolar electric field superposed in Schwarzschild background 

14 



geometry under different initial conditions. In these situations the particle starting at any 
given position and having its radial velocity inwards comes towards the central object slowly 
and as it approaches its minimum value of r (depending on r D ) and of 6 (depending upon 
the angular momentum) the particle gets bounced off in a straight line trajectory. This 
phenomenon could be a potential mechanism for the acceleration of high energy cosmic ray 
particles. In case the particle's radial velocity is directly away from the center, the particle 
spirals around the central body as depicted in Figure 3a. It is very interesting to note that 
near the object the nature of the particle trajectory remains the same (as shown in Figure 
2) whether one includes the electric field or not. However, at a very large distance from 
the central object the presence of the electric field becomes visible as the two trajectories 
differ significantly. In the case for the closed orbit the trajectory remains almost the same 
with and without the presence of the electric field as observed from Figure 3a and Figure 
3b. Thus one is lead to a very interesting conclusion that the presence of toroidal magnetic 
field totally dominates the role of the electric field in determining the motion of a charged 
particle in the presence of a strong gravitational field. 
6. CONCLUSIONS 

In order to understand the effect of curved spacetime on the electromagnetic field sur- 
rounding compact objects the trajectories of charged particles have been investigated in the 
context of accretion onto black holes. For this purpose the trajectories of a charged parti- 
cle in combined poloidal, toroidal magnetic field and rot at ion- induced unipolar electric field 
superposed in Schwarzschild background geometry have been obtained. It is found that the 
toroidal component of the magnetic field dominates the role of unipolar electric field in de- 
termining the motion of a charged particle in the presence of strong gravitational field. In 
the presence of a non-zero toroidal magnetic field an incoming particle approaches very near 
to the central object and then gets ejected away to infinity. This phenomenon could be a 
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plausible mechanism for the production of high energy cosmic ray. A rigorous magnetohydro- 
dynamical investigation including the toroidal component of the magnetic field superposed 
in curved spacetime could give rise to the production of jet like phenomenon observed in 
quasars and AGN. In the absence of toroidal field the particle is trapped into a closed orbit 
in the same way it behaves in the presence of only poloidal magnetic field except a regular 
oscillation perpendicular to the orbital motion is produced by the presence of the electric 
field. 
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Figure Captions 

Fig. la. - Projection of the trajectory of a positively charged particle (in purely poloidal 
magnetic field superposed on curved spacetime) in XY plane. The various parameters are 
E = 2.0, L = 91.77, X P = 80, Po = 2.5, (g) = -1.92, o = vr/2, (f ) = 0.3, O = 0°. In 
all figures a positive value of (^)o implies that the initial motion of the particle to be away 
from the center whereas a negative value to be towards the center. 

Fig. lb. - Projection of the trajectory of a positively charged particle in poloidal 
magnetic and unipolar electric field superposed in curved spacetime. The parameters are 
the same as used for Figure la except (^)o = —1-122 and W = 0.01. 

Fig. lc. - Same as Figure lb but without gravitational field. The parameters are the 
same except (j^) = —1-8 and L = 32.0. 

Fig. 2. - Projection of the trajectory of a charged particle on XY, XZ and YZ planes. 
The solid line represents the trajectory in the combined poloidal, toroidal magnetic and 
unipolar electric field superposed in curved spacetime while the dashed line represents that 
without the electric field. The parameters are E = 5.0, L = 17.9, Xp = 80.0, At = 4.0, p = 
6.0, O = tt/2, (g) = 0.0, W = 0.01 and O = 0.0. With the electric field (g) = -4.89, 
without the electric field (jg) = -4.93. 

Fig. 3a. - Projection of the trajectory of a charged particle in XY plane in the combined 
poloidal, toroidal magnetic and unipolar electric field superposed in curved spacetime. The 
parameters are the same as given for Figure 2 but p Q = 3.5, (^)o = +4.63, L = 41.382, and 
A T = 40. 

Fig. 3b. — Projection of the trajectory of a charged particle in XY plane in the combined 
poloidal, toroidal magnetic but without unipolar electric field. The parameters are the same 
as given for Figure 3a except (^) = +4.96. 
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